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Abstract 

Initially, we make a detailed historical survey of van der Waals forces, collecting the main ref- 
erences on the subject. Then, we review a method recently proposed by Eberlein and Zietal to 
compute the dispersion van der Waals interaction between a neutral but polarizable atom and a 
perfectly conducting surface of arbitrary shape. This method has the advantage of relating the 
quantum problem to a corresponding classical one in electrostatics so that all one needs is to 
compute an appropriate Green function. We show how the image method of electrostatics can be 
conveniently used together with the Eberlein and Zietal mehtod (when the problem admits an im- 
age solution). We then illustrate this method in a couple of simple but important cases, including 
the atom-sphere system. Particularly, in our last example, we present an original result, namely, 
the van der Waals force between an atom and a boss hat made of a grounded conducting material. 



1 



I. HISTORICAL SURVEY AND MAIN PURPOSES 



Intermolecular forces have been studied for approximately three centuries. Since 
molecules of a real gas condense into liquids and freeze into solids, it is natural to ex- 
pect that there must exist attractive intermolecular forces, a conclusion that had already 
been achieved by Newton at the end of the 17th century^. The phenomenon of capillarity - 
the ability of a liquid to climb the walls of a tube in opposition to external forces like gravity 
- was studied by the first time by Clairaut who, in 1743, suggested that this phenomenon 
could be explained if the forces between the molecules of the liquid and those of a tube 
of glass were different from the intermolecular forces between the molecules of the liquid 
themselves^. This same phenomenon was considered later on by Laplace, in 1805, and by 
Gauss, in 1830. Many other renowned scientists of that time were also involved with the 
study of determining the force law of intermolecular forces as, for instance. Maxwell and (a 
little bit later) Boltzmann. Both of them worked in the context of the kinetic theory of gases 
and based their conclusions on the available data for diffusion coefficients, specific heats and 
viscosities. Curiously, while Maxwell concluded that the intermolecular forces should be 
attractive, Boltzmann showed that repulsive forces could explain the available data as well. 
A more complete list of names involved direct or indirectly to this subject up to the 20th 
century can be found in the recent edition of the book by Israelachvili^ . 

Following a different approach, J.D. van der Waals suggested in his dissertation presented 
in Leiden in 1873^ an equation of state for real gases, given for one mol of gas by (P + 
a/V'^){V — b) = RT, where P, V and T are respectively the pressure, volume and absolute 
temperature of the gas, R is the universal constant of gases, and a and b two (experimentally) 
adjustable parameters. Parameter b was introduced by so that the finite volume of the 
molecules was taken into account: after all, the gas can not be indefinitely compressed to 
zero volume (as allowed by the equation of state for an ideal gas). On the other hand, the 
term a/V"^ is related to the existence of an attractive intermolecular force, since its presence 
leads to a smaller pressure. In other words, in a real gas the pressure is smaller than in 
an ideal gas due to the attractive intermolecular forces. These attractive forces are called, 
generically, van der Waals forces and since the seminal work of van der Waals much effort 
has been devoted to understand the nature of such forces. It is worth mentioning that J.D. 
van der Waals was awarded with the Nobel Prize of Physics in 1910. 
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In the beginning of the 20th century, M. Reinganunt^i^ described the van der Waals forces 
as a result of the interaction between the permanent electric dipoles of the molecules (he 
believed that all molecules possessed permanent dipoles). Though this is not true, since 
there are non-polar molecules, his work can be considered an important step towards the 
correct interpretation of van der Waals forces. In fact, we must distinguish three types of 
van der Waals forces: the orientation force, the induction force and the dispersion force, to 
be described below. 

Orientation forces occur between two molecules with permanent electric dipoles. Making 
a thermal average of the electrostatic interaction between two randomly oriented electric 
dipoles of moments and p2, Keesom^''' computed the van der Waals interaction energy 
between two polar molecules in a thermal bath at a given temperature and found 



where pi = |pi|, = IP2I, is the distance between the two molecules, ks is the Boltzmann 
constant and T is the absolute temperature. The minus sign in the previous equation means 
that the interaction is attractive and the subscript or stems for orientation forces. Though 
there are as many configurations which give rise to attractive forces as configurations which 
give rise to repulsive forces, Boltzmann weight [e"^^^'^'^) favours the lower energies which 
correspond to "attractive configurations". The presence of ksT in the denominator is also 
quite natural, since as T increases indefinitely all configurations (attractive or repulsive ones) 
become equally available, leading to a vanishing force. 

Induction forces are those that occur between a non-polar but polarizable molecule and 
another one which possess a permanent electric dipole (or even with a higher multipole 
moment, as for instance an electric quadrupole). Evidences that non-polar molecules indeed 
existed, together with the fact that for these substances the van der Waals parameter a was 
related to the refractive index, led Debye^"^ and others to consider this kind of forces. The 
permanent electric dipole of one molecule induces an electric dipole in the non-polar but 
polarizable one, so that we expect in this case also a behavior similar to the previous dipole- 
dipole interaction. In fact, simple arguments show that the resulting interaction energy for 
this case will be attractive and proportional to If pi is the magnitude of the dipole 

moment of the polar molecule then the magnitude of the corresponding electric field at 
the position of the second (non-polar but polarizable) molecule will be £"1(2) ~ pi/r^ and. 




(1) 
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consequently, the magnitude of the induced electric dipole acquired by the second molecule 
will be p2 = a2-Ei(2), where a2 is the static polarizability of the second molecule. Hence, 
apart from a simple numerical factor, we get for the interaction energy Uind ~ —P2Ei{2) ~ 
—a2pl/r^. Though the spatial dependence of Uindij') is the same as that for Uor{r), the 
induction force does not disappear for high temperatures, since the orientation of the dipoles 
1 and 2 are not independent as in the orientation force. Indeed, in a first approximation, the 
induced dipole 2 is parallel to the field generated by the dipole 1 at the position occupied 
by the dipole 2, which also explains the attractive character of the induction van der Waals 
interaction between a polar molecule and a non-polar one. 

The two types of van der Waals forces just described can not be used to explain the 
attraction between two non-polar molecules or two atoms, like those of noble gases (a clear 
evidence that these forces indeed exist lies in the fact that noble gases condense). In other 
words, the term a/V"^ in van der Waals equation of state will still be necessary to describe 
noble gases more accurately than the ideal gas equation of state does. However, the correct 
explanation for the forces between two atoms or two non-polar molecules, called dispersion 
van der Waals forces, had to wait for the advent of Quantum Mechanics. Due to quantum 
fluctuations, the charge and current distributions in an atom (or molecule) fluctuate and 
consequently we can think that instantaneous dipoles (or higher multipoles) exist and give 
rise to an electromagnetic interaction. These quantum fluctuations are ultimately related to 
the Heisenberg uncertainty principle, one of the most important pillars of Quantum Mechan- 
ics. In 1927, Wangi^ solved the Schrodinger equation for two Hydrogen atoms separated 
by a distance r much greater than the Bohr radius but considering the contribution of 
the instantaneous electric dipole interaction potential between the two atoms. After using 
a complicated perturbation method developed by Epsteiniiii^, he found 

C"-»M«-8.7(j^- (2) 

Three years later, Eisenschitz and London^^ and Londort^ considered the same problem 
in much more detail and used a much simpler perturbative method, refining the numerical 
factor in the previous result and relating the interaction potential for the dispersion van der 
Waals forces directly to the atomic polarizability of the Hydrogen atom. By the way, since 
the dynamical polarizability a{uj) is closely related to the permittivity e(u;), these forces 
are called dispersion van der Waals forces (name coined by London in the latter article^^). 



The expression for the dispersion interaction energy between two polarizable atoms can be 
written in the formic 

TT ( \ ^ flUJoal 
^^-^^^ = -4(W^' 
where Uq is the dominant transition frequency for the interaction and is the corresponding 

static polarizabihty of the atoms. Although the (— power law already appears in the 

first paper—, the second one, written solely by London^, is by far more cited than Eisenschitz 

and London's paper—, so that dispersion van der Waals forces are usually called London 

forces (an ordinary quantum mechanical calculation of London's result can be found in 

many textbooks, as for instance those written by Cohen- Tannoudji, Diu and Laloelfi and by 

Bransden and Joachain^^). 

For more than one decade after London's result had been established, people believed that 

the final explanation for dispersion forces was that given by London. However, experiments 

with colloidal suspensions (for a comprehensive discussion on colloids we suggest the nice 

book by J. Bergi^), made in the Phillips laboratories in the first half of the 1940's by Verwey 

and Overbeek, showed that if London's result were used there would be discrepancies between 

theoretical predictions and experimental data (see for instance the paper by Verwey''^'^ and 

the book by Verwey and Overbeek—). They noticed that in order to retrieve agreement 

between experimental data and theory the dispersion interaction energy between two atoms 

(or two non-polar but polarizable molecules) should fall for large distances more rapidly 

than Further, Overbeek conjectured that such a change in the force law was due to 

retardation effects of the electromagnetic interaction, since the velocity of light is finite. 

Retardation effects become important as the time elapsed by light to propagate from one 

atom to the other is of the order of characteristic times of the atoms, namely, 1/umn, where 

Umn are the allowed transition frequencies of the atoms. Assuming there is a dominant 

transition frequency, say Uq, retardation effects cease to be negligible for r/c > 1/ujq (in 

terms of wavelengths, this condition is written as r > Aq, where Aq is the wavelength of the 

dominant transition). Generally speaking, we can then distinguish two regimes for dispersion 

interactions, namely: the non-retarded or short distance regime and the (asymptotically) 

retarded or large distance regime. The latter is valid for r ^ Aq while the former is valid for 

ao <^ r ^ Ao, with being the Bohr radius (the condition oq <^ r is to avoid the overlapping 

of the electronic clouds of the two interacting atoms). The influence of retardation effects on 

the London-van der Waals forces was first reported by Casimir and Polder in 1946 in a very 



short paper—. Two years later, Casimir and Polder published a large paper containing all 
the details of a fourth order perturbative calculation in Quantum Electrodynamics that led 
them to their previous results^i. They showed that in the asymptotically retarded regime 
(r ^ Aq), the dispersion interaction energy between two polarizable atoms is given by 

23hc ttitts .-s 

where ai and 02 are the static polarizabilities of atoms 1 and 2, respectively. As we can see 
from the previous equation, in the retarded regime the power law of the interaction energy 
changes by a factor one, from to In this same paper, Casimir and Polder also 

showed that the dispersion interaction energy between a polarizable atom and a perfectly 
conducting plane, in the (asymptotically) retarded regime is proportional to l/z'^, where z 
is the distance from the atom to the conducting plane, instead of proportional to a 
result valid in the non-retarded regime, as showen by the first time by Lennard- Jones in 
1932^. We can understand, qualitatively, why retardation effects weaken the interaction 
between two atoms as follows. When the time taken for the electric field created by the 
fluctuating dipole of atom 1 to reach atom 2 and return to atom 1 is of the order of the 
period of the fluctuating dipole itself, the instantaneous dipole of atom 1 will have changed 
substantially from its original value so that the mutual configuration of both atoms is less 
correlated and less favourably disposed for an attractive interaction^ . It is worth mentioning 
that the orientation and induction van der Waals interactions remain non-retarded at all 
separations, only the dispersion van der Waals interaction is influenced by the retardation 
effects of the electromagnetic interactions^. The first time a transition from the retarded 
regime to the non-retarded one was observed occurred only in 1968 in the experiment made 
by Tabor and Winterton^^i^. 

In contrast to what happens to the Coulomb interaction among many point charges, 
which obeys the so called Superposition Principle, van der Waals interactions, in general, 
are not pairwise additive, as first noticed by Axilrod and Teller^^. By non-additivity of van 
der Waals interactions we mean that the interaction between two atoms is affected by the 
presence of a third one. This fact must be taken into account in the computation of the 
van der Waals force between an atom and a macroscopic body (or between two macroscopic 
bodies), since a pairwise integration using the London or Casimir and Polder forces is not 
rigorously valid anymore, except for rarefied bodies. The non-additivity of van der Waals 
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forces are ultimately related to multiscattering processes. For instance, in the case of three 
atoms, the field emanated by the first atom can reach the second one directly or after being 
scattered by the third atom. Non-additivity effects on the energy of a system may be positive 
as well as negative and are usually small (approximately less than 20%), but they can be very 
important, as for instance, in the way the atoms of rare gases are arranged in solids-. More 
details about the the non-additivity of dispersion van der Waals interaction can be found in 
Ref(s)^i2^i^ (see also Ref.'^ for a simple way of understanding this feature of dispersion van 
der Waals forces). Detailed calculations of the dispersion van der Waals interaction between 
two polarizable atoms at any separation can be found in the pedagogical paper by Holstein^^ 
and in some textbooks like those of Power'^°, Craig and Thirunamachandran^, and Salant^, 
where a careful analysis of the particular cases of the retarded and non-retarded regimes are 
presented. For an introductory discussion on the three types of van der Waals interactions 
and also the dispersion force between a polarizable atom and a conducting sphere see Taddei 
et a^. 

Though we are not going to discuss any experiment on the measurement of dispersion 
forces involving atoms and macroscopic surfaces, it is worth mentioning a few of them. 
One of the first experiments involving a beam of atoms scattering by a cylindrical surface 
was made in 1969 by Raskin and Kuscb^. In 1993, a remarkable experiment was done 
by Sukenik and collaborators^, in which for the first time the change in the power law 
between retarded and non-retarded regimes were observed directly with atoms. In 1996, 
Landragin and collaborators measured the van der Waals force in an atomic mirror based 
on evanescent waves^^. A few years later, quantum reflection was used to measure dispersion 
forces by Shimizu^J. A short but valuable description of these experiments can be found in 
the nice paper by Dalibard^. 

Dispersion forces appear not only in different areas of Physics, as in atomic and molec- 
ular physics, condensed matter physics and quantum field theory, but also in Engineering, 
Chemistry and Biology^. In Quantum Field Theory, it is closely connected to the Casimir 
effect (for a short history of this effect and its origin in experiments on colloidal chemistry 
see the introductory papers^!^ and the books— i^"— and references therein). There are even 
more bizarre situations where dispersion forces play an important role, like in the adhesion 
of geckos to the ceiling of our houses^i^ or in the formation of a thin liquid layer on ice 
that has been identified as an important element in the generation of electric potentials in 
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thunderstorms^^. A vast list of references on this subject can be found in the Resource Let- 
ter pubhshed in this journal by Milton^ and a very complete list of references on dispersion 
forces can be found in the paper by Buhmann and Welscb^. Concerning the whole history 
of intermolecular forces we recommend the excellent book by Rowlinson^°, where the reader 
can find a huge bibliography containing the relevant works written in the last three centuries 

Though the correct description of dispersion forces between two polarizable atoms sepa- 
rated by an arbitrary distance demands, somehow, the quantization of the electromagnetic 
field, which makes the general calculation very hard, the particular case of non-retarded 
dispersion forces (or simply dispersion van der Waals forces) can be computed with ordinary 
quantum mechanics avoiding completely the quantization of the electromagnetic field. In 
this article, we shall be concerned only with non-retarded dispersion forces. Particularly, we 
shall focus our attention to the van der Waals interaction energy between a polarizable atom 
and a perfectly conducting surface. Our purposes here are the following: (i) to popularize a 
simple but powerful method proposed in 2007 by Eberlein and Zietal^, and used afterward 
by them and other authors^""-^, which is extremely well suited for this kind of calculation; 
(ii) to show that the usual image method of electrostatics can be very useful in applying 
Eberlein and Zietal method (this will become clear in the explicit solutions of a couple of 
examples) and (in) to present an original result, namely, the calculation of the dispersion 
van der Waals force between an atom and a conducting "boss hat" (a conducting hemisphere 
attached to an infinite conducting plane). Our aim is to make our presentation in a level 
easily understandable to an undergraduate student with a good background in elementary 
electromagnetism and quantum mechanics. 

This paper is organized as follows. In the Section 2 we review the Eberlein and Zietal 
method showing how we can combine it with the image method. In Section 3 we illustrate 
the method by solving explicitly a couple of introductory problems, namely, the simple atom- 
plane system and the less obvius atom-sphere system. In Section 4 we treat the non-trivial 
case of an atom in the presence of a boss hat which, as far as the author's knowledge goes, 
had not appeared in the literature before. Section 5 is left for conclusions and final remarks. 
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II. EBERLEIN-ZIETAL METHOD 



In this section we review Eberlein and Zietal method^. As we shall see, this method 
has the advantage of relating the quantum problem to a corresponding classical one in 
electrostatics so that all one needs is to compute an appropriate solution of Laplace equation. 

However, with the purpose of shedding some light in the method to be explained, we 
start for convenience by making a few comments on how to compute the van der Waals 
force between two Hydrogen atoms. FIG. [T] shows two Hydrogen atoms separated by a 
distance i? = |R|, as well as other relative position vectors relevant to the problem (we have 
not included in this figure the other relative position vectors to avoid overloading it). In order 
to obtain the interaction hamiltonian, to be used in the perturbative quantum mechanical 
calculation of the interaction energy between the atoms, we write the electrostatic Coulomb 
interaction Ucoui among all the charges of the whole system. Then, we subtract from Ucoui 
the coulomb interaction between the electron and the proton of each atom and, finally, we 
make a Taylor expansion assuming that ri,r2 <^ R, where ri = |ri| and r2 = |r2|. 



— e — e 




Atom 1 ^ Atom 2 



FIG. 1: Two Hydrogen atoms whose nucleus are separated by a distance R much greater than 
their sizes. For convenience, the figure is not in scale. The relative position vectors R — ri and 
R + r2 — ri are not shown in the figure to avoid overloading it. 

The above mentioned Coulomb interaction energy is readly given by 

4vreo ln+r2i +4vreoU \^ + ^2\ |R - ri| + |R + r2 - nl / ' 
The interaction hamiltonian Hint is then easily identified as the sum of the last four terms 
on the right hand side of the previous equation. It is straightforward to show that, after a 
Taylor expansion, the dominant contribution for Hint can be written as 

rr 1 / P1-P2-3(P1-R)(P2-R) 1 

^-* = 4^\ R^ /' 

where pi = eri and p2 = er2 are the quantum mechanical electric dipole operators, and 
R = R/i?. As expected, the dominant contribution is the dipole-dipole interaction. As- 



suming the two atoms are on their fundamental states, it is easy to check that a first order 
perturbative calculation yields a vanishing result, due to spherical symmetry of the funda- 
mental wavefunctions. Performing a second order perturbative calculation one obtains^^ii^ 
an attractive interaction energy proportional to as discussed in the introduction (see 

equation (E])). 

The lesson to be learned here is the fact that in the non-retarded regime, the interaction 
hamiltonian is given by the Coulomb energy of the system, after appropriate subtractions. 
However, for more complicated situations involving continuous distributions of charges, as 
will be the case of our interest, instead of computing Ucoui by a pairwise summation of the 
coulomb interaction between all pairs of charges, we may express the electrostatic energy in 
terms of the electrostatic potential $(r), namely. 

Let us now come back to our problem, namely, that of a polarizable atom close to a 
perfectly conducting surface of arbitrary shape. This leads us naturally to the corresponding 
classical problem of a dipole near a conducting surface. The presence of the dipole induces a 
surface distribution of charges and we need to calculate the electrostatic energy of the system 
in order to determine the interaction hamiltonian to be used in the subsequent perturbative 
(quantum mechanical) calculation. As it will become evident, this task will be conveniently 
made with the aid of equation ([7]). 

The electrostatic potential $(r) satisfies the Poisson equation 

V2$(r) = , (8) 

^0 

where p(r) is the charge density, submitted to the appropriate boundary condition on the 
surface S. Assuming, initially, that we have a grounded surface, the BC imposed on the 
electrostatic potential is given by 

$(r) =0. (9) 

The electrostatic energy of the configuration is then given by equation ([7]). Solutions of equa- 
tion ([8]) can be obtained by using the Green function method^, where the Green function 
G{r,r') satisfies, by definition, the following equation 

V^G{r,r') = -6{r-r'). (10) 
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Therefore, a general solution of equation ([8]) can be written as 

$(r) = - [ G{r,r')p{r')dh' . (11) 

In order that the electrostatic potential obeys the BC written in ([9]) it suffices to impose the 
same BC to the Green function, namely. 



G(r,r': 







(12) 



res 



In terms of the Green function G{r, r') the electrostatic energy given by ([7]) takes the form 

1 



Coul 



d-^r p(r)<l>(r) 



^ jdhdh' p{r)Gir, r)p(r'). 



(13) 



One solution of equation f lTU]) is readily obtained, since this equation is nothing but the 
Poisson equation for a point charge at position r', except for a constant multiplicative 
factor. Hence, a particular solution of (fTOj) is given by 

1 



Gp(r,r') 



(14) 



47r I r — r' ^ 

However, this solution does not obey the correct BC given by f lT2|) . In order to adjust the 
BC, we add to this particular solution a solution of the homogeneous equation and write 

1 



G(r,r') 



47r|r — r' 

where ^^(r, r') satisfies Laplace equation, 

V2GH(r,r') 



0. 



(15) 



(16) 



From (fT2!) and (fTSjl . we immediately determine the BC satisfied by Gnir^r'), which reads 

1 







(17) 



res 



47r|r — r' 

Recall that all the information about the geometry of the problem is contained in (^^(r, r'). 
Let us now consider the charge density to be used in our problem. Regarding the atom as an 
electric dipole with the positive point charge at position fq and the negative one at position 
fq + h (in the appropriate moment we will take the limit h — ?■ 0), and recalling that $(r) 
vanishes on the surface, we may write 



p(r) = q 5(r - fq) - (5(r - (ro + h) 



(18) 



11 



We now substitute (fT8i) and (fT5|) into (fT3!l and after that we take the hmit h — )• 0, with 
gh — )■ d (in a moment this will be interpreted as the atomic dipole operator), we obtain 



Coul 



1 g2 1 g2 



lim , , , , X / X 

h^-o 1^ Stt^o |(ro + h) - (ro + h) 



+ 



+ 



1 



+ 



1 



ineo |ro - rol 

^2 



vreo |ro + h - rol Stt^o |ro - (ro + h)| 



+ T^[G'//(ro + h,ro + h) - G^(ro + h,ro)] 
- 7^[G'/^(ro, ro + h) - ^//(ro, ro)] \ . 



(19) 



Although the previous expression for Ucoui contains eight terms, only the last four terms 
are of interest, since only these terms contain information about the interaction between the 
dipole and the surface. Indeed, the first two terms account for the divergent self-interaction 
of the point charges at ro and tq + h and the next two terms stand for the divergent self- 
interaction of the dipole. The remaining terms can be put into a more useful form. Making 
a Taylor expansion of Gh{t^o + h, tq + h) in powers of h, it follows that 

^/^(ro + h, ro + h) - ^/^(ro + h, ro) 

= ^^(ro + h, ro) + +h ■ VGnivo + h, r') |^,^^ - G^^(ro + h, ro) 
h- V'GH(ro + h,r')L,_. (20) 



lr'=ro 



Analogously, a Taylor expansion of G//(ro, ro + h) yields 

Gniro, ro + h) - ^//(ro, tq) = h ■ VG'H(ro, r') 



(21) 



r — ro 



Hence, subtracting the irrelevant self-interaction terms mentioned before and using equations 
( 120|) and ( 12T]) . the relevant interaction hamiltonian operator to be used in the perturbative 
quantum mechanical calculation of the van der Waals force between the atom and the 
(grounded) conducting surface is given by 



Hi^t = lim-^(d-V') G^(ro + h,r')-G'H(ro,r') 



r'=ro 



— (d-V')(d-V)G^(r,r') 



(22) 

r=r'=ro 

It is worth emphasizing that the atomic dipole moment that appears in this expression is a 
quantum operator. In first order of perturbation theory the desired non-retarded interaction 
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energy between the atom and the conducting surface, denoted by Unr, is just the quantum 
expectation value of the above expression, namely. 



1 

f^iVi?,(ro) = {Hint) = 7^ (rfmC^n)VmV^GH(r,r') 
ZEq — 



(23) 

m,n=l r=r'=ro 

where we used the fact that the only operators in the above expression are and dn 
{Gh{i^, r') is a c- number). For the sake of simplicity we shall always work with an orthonor- 
mal basis, for which we can write 

(dmdn) = 5mn{dl^) ■ (24) 

Combining the two last equations we finally obtain 

1 3 

Unr{vo) = — $^(rf^)5^5;GH(r,r')|r . , (25) 

m=l 

which is precisely the expression obtained by Eberlein and Zietal^^. This method has the 
advantage of relating the quantum problem to a corresponding classical one in electrostatics. 
Its remarkable simplicity consists in the fact that to obtain the non-retarded van der Waals 
interaction energy of an atom near any conducting surface one must find only the homoge- 
neous solution of Laplace equation, Gh, corresponding to that geometry. In other works, 
one must solve the classical problem defined by equations (fT6|) and (!T7|) . These equations 
are, except for constants, precisely those that yield the electrostatic potential of the image 
charges for the problem of a charge at position r' in the presence of the surface S, (if the 
problem admits a treatment by the image method). 

Indeed, denoting by 0(r) the electrostatic potential of that configuration, we may decom- 
pose it into the sum of the potential of a charge plus the potential of the image charges, 
denoted by 4>i{r). Since the potential of the single charge satisfies Poisson equation of the 
problem, 0j(r) obeys Laplace equation. Together with the boundary condition, the equations 
satisfied by 0.j(r) are the following 

VVi(r) = (26) 

r H+<^^(r)l =0. (27) 

From a direct comparison between equations fl26|) and fl271) and equations (|T6l) and f fT2|) . it 
is straightforward to make the identification 

GH(r,rO = ^. (28) 
13 



Note that the dependence of the rhs of previous equation on r' is imphcit, since the image 
charge depends on the position of the physical charge - recall that the physical charge is 
located at r'. 

In conclusion, the image method is a very useful tool in order to find the homogeneous 
solution Gh which, in turn, is the only function needed to perform Eberlein-Zietal calcu- 
lation and obtain the quantum non-retarded dispersive interaction between an atom and a 
conducting surface S of an arbitrary shape. In the next sections we will apply this procedure 
for different geometries. 

III. INTRODUCTORY EXAMPLES 

In this section, we apply Eberlein and Zietal's method just discussed in introductory 
examples, namely, for an atom close to a grounded conducting plane and an atom close to 
a grounded conducting sphere, where the image method can be employed with no difficulty. 
In the latter case, it is possible to solve for a non-grounded isolated sphere as well. However, 
in this case, appropriate modifications of the method outlined in the previous section are 
needed, since (^//(r, r') no longer satisfies the boundary condition (|T2l) . 

A. Atom close to a grounded conducting plane 

Consider a polarizable atom at position vq in the presence of an infinite conducting plane 
located at 2; = 0. As outlined in last section, all we need to obtain the dispersion van der 
Waals interaction energy for such a system is to find out the function Gj:/(r, r') associated 
to it. FIG. [2] shows the charge q at position r' = {x',y',z') and the conducting plane, as 
well as the image charge qi which, in this case, is simply given by qi = —q and is located at 
position r/ = {x' ,y' , —z'). 

The electrostatic potential at r created by the image charge —q located at position r/ is 




(29) 



where r-=r' — 2z'z. Hence, from equation f l28|) . we readly obtain 




(30) 
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FIG. 2: Point charge q near an infinite conducting plane and its image. 

Now we are ready to use Eberlein and Zietal's method. Substituting the previous expres- 
sion for Gnir, r') into equation (p5|) . we have 

3 



1 f 

m=l ^ 



1 



r=r '=ro 



(31) 



[{x - x'y + {y- y'f + {z + z'YY/-^ 

The derivatives in the previous expression can be easily computed. For instance, for the 
first coordinate, m = x, we get 



[{x - x'f + {y - y'Y + {z + z'YYl^ 



1 



r=r '=ro 



(32) 



An identical result is valid for m = y, while for m = z, the following result is obtained 



4ko| 



(33) 



[{x - x'y + {y- y'Y + {z + z'fyl^ 
Substituting the previous results into equation fl3T|) . we finally obtain 

Ua,{zo) = - '\ ' . 34 

This is the well-known interaction between an atom and an infinite conducting plane in the 
non-retarded regime firstly obtained in 1932 by Lennard-Jones-- (see also^"^). 



B. Atom close to a grounded sphere 

Now let us consider an atom in the presence of conducting grounded sphere of radius R 
and center C. The corresponding electrostatic problem we need to solve is that of a point 
charge q at position r' = {x', y', z') in the presence of the conducting sphere. The image 
method for this problem tells us (see, for instance, Griffith's textbook—) that we have to 
put an image charge qi = —fjq at position r/ = where r' = |r'|, as sketched in FIG. El 
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FIG. 3: Point charge q near a conducting sphere of radius R and its image Qi. 



Hence, the potential of created by the image charge at a generic point r {r > R) is given by 

'^^^'"^ = 47r£o|r-r/| = ~4w^^-r/| " ^^^^ 

Substituting this expression into equation (128|) we get the homogeneous solution Gnir^r') 
for this configuration, 



GH(r,r') 



R 



n ■ (36) 

q 47rr'|r — | 

Hence, all we have to do is to apply Eberlein-Zietal formula, i.e., equation (1251) . For sim- 
plicity, we shall consider an isotropic atom, a hypothesis which allows us to write 



(37) 



For convenience, we can orient the axis so that the atom is located at fq = {0,0, zq). Note 

that with this assumption there is no loss of generalization, since the problem exhibits 

spherical symmetry. The calculation is lengthier than the previous one for the atom-plane 

configuration but it still involves only elementary derivatives, so that we shall show only the 

main steps, leaving for the readers the verification of the intermediate ones. 

-1 



Let us compute d^dx I r'|r — r/ 



X + y + z 



r=r'=(0,0,zo) 
1/2 



where 



/ y ( R^ ,y / 
\^-—2^) +\y-—2y) 



n 1/2 



With this purpose, note initially that 



X 



ii 

.12 



{x-^,x'y+{y-^,y'y+{z-fr-2z') 



i?2 



3/2 



(38) 



(39) 
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By computing the derivative of the previous expression and after that evaluating the 
result at r = r' = (0, 0, Zo), it is straightforward to show that 

i?2 



did: 



(40) 



r=r '=(0,0,^0) ^0 (1 " ^o)^ 

An identical result is obtained for the coordinate m = y. In the case of coordinate m = z 
it can be shown with an analogous but lengthier calculation that 



M J« J» ^ 



2^' 1 

+ ^.,0.3 ■ (41) 

r=r'=(0,0,^()) 



zi (1 - nyzif zt (1 - w/zif ■ 



Collecting all the previous results and substituting them into (125|1 . we finally obtain the 
dispersion van der Waals interaction energy between an atom and a grounded conducting 
sphere, 

f/aps(^0, R) - "24^^^ I ^6 _ ^,/^,^3 + ^4 (1 _ ^2/^2)2 ] ' (42) 

For later convenience, we rewrite the previous result in terms of R and the distance between 
the atom and the sphere, zq — R, which we shall denote by a. Substituting a = zo — R in 
(142|) we get, after a few rearrangements, 

TT I (d^) / 4 alR \ . 

^"^^^^"'^^ - ~2A^ [ (2 + a/Rf + (2TWJ ' ^ ^ 

For an atom with a dominant transition frequency, we may still cast our result in terms of 
the static polarizability of the atom. Recalling that the static polarizability for an atom in 
its ground state is given by— 

2^^ (44) 

where ojno is the transition frequency between the n-th state and the ground state, and d„o 
is the corresponding transition dipole moment, for an atom with a dominat transition, say 
between the fundamental state and the first excited one, we have 

Since for an atom with this dominant transition (d^) = |diop, equation (H3l) reduces to 

TT I T-.\ /iwioct f 4 alR 1 

^-^"'^^ = -16^ 1 (2 + alRf + (2TWJ • ^ ^ 

Some comments are in order: (i) the previous result was obtained by the first time by 
Taddei and collaborators^. However, the agreement of our result, given by equation ()i6|) . 
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and the result obtained by these authors is up to a numerical factor of 3. A discrepancy of 
a numerical factor between both results was expected, since in Ref.'^ the authors employed 
the (semiclassical) fluctuating-dipoles method, which is not expected to provide the correct 
numerical factors, though it gives the correct behavior of the interaction; (ii) the atom- 
sphere system had been discussed before by many authors^^"— for spheres with different 
properties and for regimes other than the non-retarded one and recently has been a subject 
of great interest^"—. In fact, the result expressed in equation fH3l) was reobtained by 
Buhmann as a particular case of a more general discussion^; (in) equation P6l) is valid for 
any values of R and a, provided the conditions for the non-retarded regime remain valid. It 
is easy to show that, in the limit R — > oo, with finite a, equation fl43p reproduces the non- 
retarded interaction energy for the atom-plane system, namely, Uags{zo) — )■ — (d^)/[487r£:o a^], 
in agreement with equation flM|) if there we write {d"^) + {dy) + 2{d'j.) = (4/3) (d^). 



C. Atom close to an isolated conducting sphere 

Let us consider now with a neutral perfectly conducting isolated sphere. This case differs 
from that of an atom close to a grounded sphere because now Guij, r') does not satisfy the 
boundary condition written in ( |T2|) . Since the sphere is not grounded anymore, the presence 
of a dipole changes its potential. Hence, the first thing we have to do is to find out the 
electrostatic potential on the surface or the sphere. With the aid of the image method, it is 
possible to show that the potential induced on an isolated sphere by a dipole at position fq 
is given by^ 

Therefore, the BC to be satisfied by the electrostatic potential in this problem is 

d ■ rp 

47r£:o|ro|^ 

Substituting the above result into equation (fTTIl . we obtain 



*«||r|=H=7Zr^- (48) 



G{R f , r')p(r') dh' = -^-^ . (49) 



47r ro 



Using the same charge density as before, namely, that given by (1181) . and following the same 
procedure employed to go from equation (|T9|) to equation fl22l) . we get 

d-V'G'(i?r,r')| , =i^r^. (50) 

^ ^lr'=ro 47r|ro|3 ^ ^ 
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Hence, in order to the boundary condition ( HHj) be fulfilled, we may impose the following 
BC to the Green function, 

VG{t, r' 



= TT^s ■ (51) 



r ' = ri 







Apart from this condition, the Green function must also obey Poisson equation (iTOl) . As it 
will become evident in a moment, a convenient way to obtain the desired Green function is 
write it as the sum of the Green function of the previous case (with the grounded sphere), 
with an extra term G^^\r, r') to be determined, 

^(r, r') = - -— ? + G« (r , r') . (52) 

47r|r — r I 47r|r | |r — r^'l 

The first two terms on the right hand side of the previous equation, together, satisfy the 
Poisson equation in the presence of a grounded sphere, so that these two terms (when added), 
vanish on the surface of the conductor. Therefore, in order to solve the problem with an 
isolated sphere, we just impose that G^^^{r,r') must satisfy the Laplace equation as well as 
the boundary condition flSTj) . to wit. 



V^G'^'\r,r') = (53) 
V'G(^)(r,r' 



= T?^3 ■ (54) 



••0 



Except for a constant factor, the right hand side of fl54j) is identified with the electric field 
created by a point charge at the origin. Since the electromagnetic is minus the gradient of 
the electrostatic potential, G^^^Rr, r') is naturally identified with the electrostatic potential 
created at a point r' by a point charge at the origin. Note that this is compatible with 
equation since this equation must be satisfied only in the physical region of the problem 
at hand, namely, the region outside the sphere. Hence, we may write 

where /(r) must assume the unit value for points belonging to the surface of the sphere, in 
order to satisfy equation (15^ . Using the fact that the Green function must be symmetric^ 
by the exchange r < — > r', and since /(r) = 1 for |r| = R, we are led to the result 

Ga)(r,rO = , , . (56) 
^ ^ 47r r r' ^ ^ 
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Inserting ( 156|) into (152|) we see that the G_f/(r, r') function for the present case is given by 

1 



GH{r,r') = G{r,r': 



47r|r — r'l 
R R 

77 + 



(57) 



47r|r'| |r — r/| 47r|r| |r'| 

where r/ = ^r', as in the grounded case. 

As we have seen, some care must be taken with problems involving isolated conductors, 
because the boundary condition will depend on the charge distribution of the system. We 
should also emphasize the convenience of having used the solution of the problem with a 
grounded sphere as a step in the search of a solution to the problem with an isolated sphere. 

As one may readily verify, equations f|T3|) -f l25|) do not depend on the boundary conditions 
satisfied by G'//(r, r') and, as a consequence, Eberlein-Zietal expression ( |25|) is still valid 
in this case. Using the previous results for the grounded sphere, as well as the following 
relations involving G^^\r,r'), 



^^^^ ; 4vr|l'i 



r=r'=(0,0,2()) 



R 



47r|r| Ir' 



dydl 



R 

4^ 



R 



47r|r| Ir' 







r=r'=(0,0,zo) 



(58) 



r=r '=(0,0,^0) 

we finally obtain the dispersion van der Waals interaction energy between an atom and an 
isolated conducting sphere. 



Uais{zQ, R) 



{(p) r4i?3 



+ 



R 



R 



247reo I zl (1 - R^ zlf (1 - / z^f 



(59) 



As in the previous section, we can also write last expression in terms of R and the distance 
from the atom to the surface of the sphere, zq — i?, denoted by a. Doing this, we have 

(d^) [4 , a/R a^R^ 



Uais{a,R) 



+ 



(60) 



247reoa3 {{2 + a/Rf {2 + a/Rf {I + a/Rf 
in agreement, up to a numerical factor of 3, with the result obtained by Taddei and 
collaborators^^. The second term on the right hand side of last equation is, in absolute value, 
greater than the third. Therefore, the interaction of an atom with an isolated conducting 
sphere is always attractive. Since the only difference between the grounded and isolated 
cases is the last term present in equation (1^ . we conclude that the attraction is stronger 
in the case of a grounded sphere. This can be physically understood as a consequence of the 
charge acquired by the grounded sphere. 
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We finish this section by taken the interesting hmit i? — )■ 0, but with AneoR^ — )■ where 
as is the (finite) polarizability of a very small conducting sphere. In this limit, the previous 
equation reduces to 

lim Uais(0',R) = ^ — ^ = — ■, (61) 

where in the last step we assumed that the transition from the fundamental state to the first 
excited one is dominant. Note that the result is a London like dipole-dipole interaction, as 
expected. 



IV. ATOM CLOSE TO A CONDUCTING BOSS HAT SURFACE 

Having solved the simple cases of the last section, we are now in position to solve a more 
interesting case, namely, an atom near a the conducting surface with the shape of a "boss 
hat" . This surface consists of a conducting spherical hemisphere with radius R together with 
an infinite conducting plane. The geometry in question and the necessary image charges to 
the problem are sketched in FIG. HI 

{pii^ziiP' ^ 4>\ p'-I^z'^ 

(^T^^rp^p', 0', —^72^777^0 

Qi3 = -q {p',^',-z') 

FIG. 4: A physical point charge q close to the conducting boss hat surface and the three necessary 
image charges located in the non-physical region. 

To begin with, consider a charge q at position r' in the presence of the conducting boss 
hat. Now put a charge qn = —frQ at position r/^ = ^r'. In cylindrical coordinates these 
equations take the form qn = --j=M=q and r/^ = (jt^jtiP' , (j)' , ji^jnz'), since r' = |r'| = 
-)- z'"^. As used in the atom-sphere case, this pair of charges, q and qn, furnishes a null 
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potential at the spherical part of the conducting surface. But the potential generated by 
these two charges does not yet satisfy the desired BC at the plane part of the conductor. 
Therefore, we must introduce two more image charges, one, with charge = —qn, being 
the mirror image of qn and the other with charge g^s = — g, being the mirror image of the 
physical charge q. The addition of these two charges, and g^a, leads to a null potential 
at the plane z = but with the advantage of not disturbing the null potential at the 
hemisphere, since by symmetry qi2 is precisely the image charge of g^s with respect to the 
sphere. Consequently, the four charges, the real charge q plus the three image charges qn, qi2 
and qi3 lead to an electrostatic potential which is zero on the boss hat conducting surface. 
The positions of the four charges, namely, r' = {p',Lp',z'), r/^ = {-^j^jnP' , (p' , J70j72z'), 
r/2 = i p'if^'i P', 0', - p'if^ri z') and r,^ = {p', ip' , -z') are indicated in FIG. H 

With this image configuration, the potential generated by the image charges is just a 
superposition of the potentials created by the charges g^i, qi2 and 9,3, namely, 

M^) = J—- \ . > I + . , I + . > I \ ■ (62) 

Using the previous expressions for qn, qi2 and gjs, as well as for r/^, r/2 and r/3, and defining 
the quantities ^(r,r'), ,^^(r, r') and ^4.(r, r') as 

C_(r,r') = |r-r;j ; ^+{r,r') = \r - r-^\ ; ^(r, r') = |r - r/g] , (63) 

the function Gni^, r') for the boss hat case can be written in the form 



where 



^(r,r') = + p2 + {^z' + zf - 2p'p cos(0' - 0) , (65) 

^±(r,r') = V^V'^ + (P'^ + z'^yp^ + [(P'^ + ± R'^zf - 2i?2(p'2 + z''^) p' p cos{(j)' - 0) (66) 

and we also used equation (128|) . 

The dispersion interaction energy of an atom at a generic position (po, (po,zo) and a boss 
hat conducting surface, UabhiPo, <Po: Zq), can then be obtained from equation (1231) . After a 
lengthy but straightforward calculation, it can be shown that 

Uabh{po,zo) = -TT- — 3\{dl)'^p{po,Zo) + {dl)E^{po, zo) + {dl)E^{po,zo) \ , (67) 



GAireoZQ 
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where 

Sp(Po, ^o) 



1 - SRz^ 



[{pl + zl + R^f - AR^plf/^ 
1 1 



{pl + zl - R^y- 



with 



2 + 



° 1 {pl + zl - R^f [{pl + zl + R^f - 4i?2p2]3/2 

C{R,Pq,Zq) 



mzl 



{pl + zl - R'Y 



R' + z' + 



[{pl + zl + R^y - 4i?Vg]^/^ 



C{R,pQ,zo) 



R'pI -10p^4 - lOptR'zl - 10i?Vo + 8plR'zl - z', + 2p<',zl+ 



+ 8plz^ - SQplRh^ + 10p2i?6 
- 5plz',{zl + pl)\izl + plf - 



2 2 

Po^o 



(71) 



As expected, the interaction energy does not depend on (p, due to the axial symmetry of the 
system. Further, one may immediately recover the atom-plane result, given by fl34p . just 
taking i? = in the previous equations. In FIG. E] we plot the interaction energy given by 
(1671) multiplied by R^ (apart from a constant factor) in terms of z^/R for the particular case 
where the atom is on the OZ axis (po = 0) and its atomic polarizability in this direction is 
dominant, i.e., {dl) ^ {dV), {dV). Note that only the function is necessary. 




FIG. 5: Dispersion interaction energy between an atom and a conducting boss hat surface as a 
function of zq/R with the atom on the OZ axis and assuming that {d^) ^ (f^p), {d'^)- The graph 
is plotted in arbitrary units. 



23 



It is interesting to analyze the curvature effects on the interaction between the atom and 
the boss hat surface by comparing the interaction for this case, given by fl67|) with that for 
the atom-grounded sphere case, given by fl42p . To be consistent, we shall now consider an 
isotropic atom in equation flB7|) . since we made this assumption in obtaining fH21) . However, 
we shall compare these two expressions only up to third order in {zq — R)/R. Making, then, 
a Taylor expansion of equations f H2|) and ( 167|1 and maintaining only terms up to third order, 
we obtain for the respective expressions of Uags{zo,R) and Uabh{zo, R): 



Uagsi^o, R) 



and 



Uabh{zo, R) 



1 



zo-R ^ (zo-Rf 7izo-Rf 



R R^ 



zo-R^izo- Rf 



7(^0 - Rf 



+ 



+ 



(72) 



(73) 



R ' R^ 8i?3 
Comparing last equations we see that they coincide up to second order in {zq — R)/R (up 
to order (2:0 — RY /R^, the interaction of an atom with a boss hat surface is the same as 
that of an atom with a sphere). This is reasonable, since an atom very close to a boss hat 
surface can not distinguish it from a sphere. As the distance between the atom and the boss 
hat increases, the differences between the two surfaces become apparent. The comparison 
between Uags{zo, R) and Uabh{zo, R) expanded up to third order in (2:0 — R)/R and the exact 
result for the interaction in the atom-grounded sphere case is illustrated in FIG. |6l 




6 8 10 Zo/R 



FIG. 6: (Color online) Exact expression for Uags (blue solid line), given by (j42p and expansions up 
to third order for Uabh (dashed red line) and Uags (dotted-dashed black line), given respectively by 
r2]) and ([73|) . as functions of zq/R. Both axis are in log scale and we are using arbitrary units. 
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V. CONCLUSIONS AND FINAL REMARKS 



We discussed here a method for computing the van der Waals interaction between an atom 
and a conducting surface of arbitrary shape introduced recently by Eberlein and Zietal^ 
and apphed it in a variety of situations, from introductory examples, as the cases of an 
atom interacting with an infinitely grounded conducting plane, or with a conducting sphere 
(grounded or isolated), to the more complex case of an atom interacting with a conducting 
boss hat surface. This method has the advantage of quickly leading us to a correspond- 
ing classical problem in electrostatics, which allows us to use well establihed methods for 
solving Laplace or Poisson equations as, for instance, those based on Green functions or on 
the image method. Particularly, we have shown explicitly how the image method can be 
appropriately employed jointly with Eberlein and Zietal's method (when the problem admits 
an image solution). We should emphasize, however, that the use of the image method is 
not mandatory. In fact, in the pioneering work on this method^, the authors discussed the 
problems of an atom interacting with an infinite conducting semi-plane, and an atom inter- 
acting with a conducting cylinder using the Green function method. Also, in further papers, 
this method was employed without the use of the image method, as in the calculation of the 
non-retarded interaction between an atom and a dielectric slab^^ and in the computation of 
the non-retarded interaction of an atom and an infinitely conducting plane with a circular 
hole^ (this same system as well as the atom-disk system were discussed with the aid of 
the more involved Sommerfeld's image methoct^ in a recent article^). Hence, it is not an 
abuse to say that Eberlein and Zietal's method is due to its simplicity and at the same time 
its power, with a very favorable cost-benefit analysis for undergraduate and also graduate 
students that are beginning to study non-retarded dispersion forces between atoms and con- 
ducting surfaces. A lot of other systems can be handled by this method. We challenge the 
interested reader, for instance, to reobtain the non-retarded force between an atom and a 
conducting wedge with aperture angle equal to vr/ra, with n a positive integer, first obtained 
by Mendes and collaborators^^ (the retarded Casimir-Polder interaction between an atom 
and a conducting wedge had already been calculated by Brevik and collaborators^). 

There are many other methods of computing dispersion forces between atoms and macro- 
scopic bodies which are much more general than the one discussed here, in the sense that 
they consider all the distance regimes, thermal effects as well as all kinds of materials, not 
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only perfectly conducting ones. It is worth mentioning, for instance, the famous results 
obtained by Lifshitz in 1956''^, and generalized by Dzyaloshinskii, Lifshitz and Pitaevskii^ 
a couple of years later, that since then have been applied with an enormous success. Their 
results predicted, for instance, the variation of the thickness of thin superfluid helium films 
in a very good agreement with the experiments^. However, the discussion of more general 
methods of computing dispersion forces would lead us very far from the main purposes of 
this pedagogical article. 

We finish this article by mentioning that dispersive forces are still the subject of an 
intense research, mainly in connection with repulsive dispersive forces, which would be 
very important in a variety of situationa^i^i^"— , to mention just a few recent results on 
repulsive forces. We should mention that Feinberg and Sucher— had already shown that 
an electrically polarizable atom and a magnetically polarizable one repelled each other. This 
result has been reotained in a simple way in Ref(s}^'^ (in connection with these results, 
see also the papers by Boyer— ). It is worth mentioning that the interaction between an 
electrically polarizable atom and a magnetically polarizable one has been discussed in a 
more general context, namely, with the atoms embedded in a magneto-dielectric medium by 
the authorsS'^. The dispersion interaction between a ground state atom and a corrugated 
surface was firstly discussed by Messina et a/^°, where the so called scattering approach 
was employed. The existence of a lateral Casimir and Polder force gave rise to several 
interesting proposals of experiments on dispersive forces, some of them including Bose- 
Einstein condensates near periodic gratings. We hope this paper can motivate the readers 
to turn their attention to such an interesting and interdisciplinary subject as the dispersion 
forces are. 
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